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Показано влияние продолжительности нагрева и параметров тепловыделения на максимальную температуру активного элемен- 
та при отсутствии отвода теплоты. Точность расчета температурного поля в активном элементе зависит от невязки дифференци- 
ального уравнения теплопроводности и числа Фурье. Установлен диапазон входных параметров, при которых сложная двумер- 
ная задача теплообмена сводится к одномерной. 


В [1] для обоснования точности инженерного 
расчета нестационарного температурного поля в 
активном элементе конечных размеров было пред- 
ложено ввести в практику критерии качества расче- 
та. Они позволяют провести проверку результатов 
аналитического расчета на раннем этапе моделиро- 
вания теплового процесса в конкретном элементе. 
Проверка состоит в подстановке расчетных значе- 
ний в исходные дифференциальные уравнения и 
краевые условия исследуемой задачи. После этого 
этапа рекомендуется провести сравнение с опытны- 
ми или другими надежными данными и приступить 
к самому процессу моделирования в широком диа- 
пазоне изменения параметров исходной задачи. 

Целью данной работы является выяснение вли- 
яния входных параметров на тепловое состояние 
нажимной плиты турбогенератора. В таком актив- 
ном элементе распределение удельных тепловых 
потерь подчиняется следующей зависимости 
Ро(Х, У, Ро) = Ро 0 • Щ (X) • ]У 2 (У) • ехр(-5Ро), ( 1 ) 


где ЩХ) = ехр(-АХ), ЩУ) =1 +МУ + ИУ 2 . (2) 

Ѳ(Х У Ро) = У У Т ' ІНі ’ V™ ’ Ро) ^і (У ,„ , У) (3) 
К п (р п )К 22 (у ш ) 


Здесь ц„, у т - собственные числа. Они находят- 
ся из трансцендентных уравнений 


сі§ц = 


р 1 -Ві,Ві 2 
ц(Ві, +Ві 2 ) ’ 


( 4 ) 


сі %уК -- 


У 1 ~ Ві,Ві 4 
7(Ві, + Ві 4 ) 


( 5 ) 


К,(р п ,Х) = р п С08 Р П Х + ВІ 2 ЗІП д X, 
Х 2 (.У т >У) = У т С08 Ут У + Ві 4 8ІП у т У; 

|д 2 + Ві 2 2 +(д 2 -ВІ 2 2 )х 1 
^п(Ц„)4| х 8т^д +в . 2(1 _ соз2л) |, 
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т, ( Д, , У„, , Ро) = Ро 0 Рфр„)Р 2 (у т , К)Р,(р„ , У т , Ро), 
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УУ У • 

ЩЯ) = 1 + МЯ + ОЯ 2 ; 

К 2 (Ут > Я) = у т соз у т Я + Ві 4 зіп у т Я; 


Р,(Р»>УтХо) - 
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г 2 -|ехр(-5Ро) -ехр[-(д 2 +у т 2 )Ро]|. 

Нп + Ут~$ 

При отсутствии охлаждения Ві, 2 _ 3 4 — >0 собствен- 
ные числа находятся из уравнений вида: 
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8Іп /л - 0; уЦ, = 0, /л 2 = п,...,ц п = (я-1)7г; 
зіпуі? =0; у 1 = 0 ,у 2 = л/ К,..., у т =(т-Ѵ)п/ К. 

Именно в этом случае наблюдается наибольшая 
погрешность расчетов [1]. 

Исследование температурных полей в нажимной 
плите турбогенератора при распределении удельных 
тепловых потерь, согласно (1), (2) и отсутствии ох- 
лаждения показало, что изменением температуры в 
направлении оси У можно пренебречь. При этом ока- 
зывается, что на точность расчета максимальной тем- 
пературы в активном элементе оказывает выбор чис- 
ла членов кнр двойной суммы ряда (3) (табл. 1,2). 


Таблица 1. Влияние времени, тепловыделения и числа чле- 
нов ряда (3) на максимальную температуру в на- 
жимной плите турбогенератора 


Ро 

Ро 

і 

І 

Ѳ 


9Ѳ/д Ро 

6 2 Ѳ/дХ 2 

д 2 Ѳ/дУ 2 

В(Ѵ>10- 3 , 

0,01 

137,23 

26 

3 

1,218 

0,0002 

112,34 

-24,82 

-0,068 

-0,318 



13 

13 

1,200 

2,0459 

110,70 

-24,48 

-0,005 

-1,357 



4 

29 

1,201 

1,9446 

110,77 

-24,49 

-0,017 

96,137 

0,05 

126,68 

26 

3 

4,991 

0,0002 

81,15 

-45,25 

-0,276 

-0,293 

0,10 

114,62 

26 

3 

8,513 

0,0002 

61,44 

-52,72 

-0,463 

-0,266 


Примечания: 

1) Исходные данные: к=30, р=30; Ро 0 =112 ВІ П ц=0,001; В=7,5; 
5=Ы=2,0при М=Щ 0=-7/й 2 ; р,=у=0,00І;Х=0,0; У=0,5Р. 

2) 1;=дѲ/дРо—(д 2 Ѳ/дХ 2 +д 2 Ѳ/дѴ 2 +Ро(Х,Ѵ,Ро)) —невязка в диф- 
ференциальном уравнении теплопроводности. 

3) ва,])=1/(р 2 +у 2 -5){-5ехр(-5Ро)+(р 2 +у 2 )ехр[-(р- 2 +у^ 2 )Ро]} 


Таблица 2. Влияние условий охлаждения на температуру в на- 
жимной плите турбогенератора (в точке Х=0,01; 
Ѵ=0,0Ш) при Ро 0 =112 ; Р=7,5 ; М=0, 0=-Щ, 
5=Ы=2,0; Ро=0, 1; Ро(Х, Ѵ,Ро)=89,87-расчетпо (1), (2) 


Вариант 
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Р 

в 1 1,23,4 

Д\=У\К 

п 

т 

Ѳпо (3) 

4 

1 

10 

10 

0,00001 

0,0001 

1 

1 

5,543 

39,8 

2 



0,0001 

0,0001 

1 

1 

-8,021 

162,3 

3 



0,0001 

0,001 

3 

2 

7,162 

0,2587 






1 

3 

6,765 

-0,0839 

4 



0,001 

0,001 

1 

3 

6,775 

0,0239 

5 

20 

20 

0,001 

0,001 

3 

2 

7,151 

-0,7876 






5 

6 

6,748 

-0,0056* 


Примечание: * — і=5, ]=6 


При изменении параметров Ро<0,05; 
0,001<Ві 1 2 3 4 <0,01; 0,075<Л<7,5; 0,2<5<2,0; 

0,2</Ѵ<2,0 при М=\/К, 0 = — 1 /IV рекомендуется в 
инженерной практике использовать более простую 
аналитическую зависимость 


и(Х,У,Ро) 


^ Т х {р п ,У,Уо)Кфц п ,Х) 

к *„оо 


( 6 ) 


Г(ц„ХРо) = Ро 0 Д(д,)Щ 2 (У)С( Н ,,Го), 


Р(п„, Ро) 



— [ехр(-5Ро) -ехр(-д 2 Ро)]. 
о) 


На рисунке показано изменение температуры в 
направлении оси X, из которого видно, что макси- 
мальное отклонение (до 5 %) температур, вычис- 
ленных по (3) и (6), отличаются на поверхности ак- 
тивного элемента. Во внутренних точках нажим- 
ной плиты это различие температур отсутствует. 
Расчет проведен на основе входных параметров: 
симметричные условия охлаждения Ві=0,01; 
Ро=0,001; Л=0,075; М=\/К, 2)=-1/7С; Ро 0 =1 12, 
У=0,5К- $=N=2,0. 



Рисунок. Расчет температуры в нажимной плите турбогенера- 
тора: 1) поур. (3); 2) поур. (6) 


Таким образом, при тепловом расчете активно- 
го элемента электрической машины рекомендуется 
учитывать влияние входных параметров (геометри- 
ческие размеры, электротеплофизические свой- 
ства, условия охлаждения, распределение внутрен- 
них источников теплоты) на точность определения 
максимальной температуры. Она существенно за- 
висит от качества аналитического расчета, опреде- 
ляемого значениями невязки дифференциального 
уравнения теплопроводности и числа Фурье. 


Сравнение вариантов 1 и 2 табл. 2 показывает, 
что не учет влияния входных параметров может 
привести к неправильному определению макси- 
мальной температуры (вариант 2, 0=-8,О21). Уста- 
новлено наличие сочетания числа слагаемых (/,/), 
при которых наблюдается минимальное значение 
невязки | по дифференциальному уравнению теп- 
лопроводности и экспоненциального множителя 
В{іф). В варианте 5 при /=5, /= 6 (см. табл. 2) получе- 
но точное значение температуры 0=6,748, соответ- 
ствующее малой величине невязки относительно 
вариантов 3^1. 


Обозначения 

Ѳ(Х,У,Уо)=(Т(х,у,т)-Т ж )/ Т м - безразмерная тем- 
пература; Т(х,у,т), Т ж , Т м - соответствующие темпе- 
ратуры, К; Ро(Х,У,Ро)=^х,у ) г)-й 2 /(Я 1 Г„,) - число По- 
меранцева; /V, Д М, В - коэффициенты функций 
распределения, учитывающие неравномерность 
внутренних ис точников тепла; Ві 12 =аі 2 'й/Я„ 
Ві 3 4=а:3, 4-й/лЙД^ -число Био; Х=х/Ъ, У=(Х,/Х у )' п у/Ъ, 
К=^Х^' П Н/Ъ - безразмерные координаты; 
Ъ, Н - геометрические размеры, м; X^X У - коэффи- 
циенты теплопроводности, Вт/(м-К); а, 234 - коэф- 
фициент теплообмена на соответствующей охлажда- 
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емой поверхности, Вт/(м 2 -К); Ро=аг/й 2 - число 
Фурье; к, р - принятое ограниченное число членов 
каждого ряда при расчете конкретного варианта рас- 
сматриваемой задачи; і, } -число членов ряда при ко- 


торых наблюдается минимальная невязка в уравне- 
нии теплопроводности. Индексы: ж - охлаждающая 
среда; м - масштаб. 
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